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We examine harmonic oscillator defects coupled to a photon held in the environs of an 
optical hber. Using techniques borrowed or extended from the theory of two dimensional 
quantum helds with boundaries and defects, we are able to compute exactly a number 
of interesting quantities. We calculate the scattering S-matrices (i.e. the rehection and 
transmission amplitudes) of the photons off a single defect. We determine using techniques 
derived from thermodynamic Bethe ansatz (TBA) the thermodynamic potentials of the 
interacting photon-defect system. And we compute several correlators of physical interest. 
We hnd the photon occupancy at hnite temperature, the spontaneous emission spectrum 
from the decay of an excited state, and the correlation functions of the defect degrees of 
freedom. In an extension of the single defect theory, we hnd the photonic band structure 
that arises from a periodic array of harmonic oscillators. In another extension, we examine 
a continuous array of defects and exactly derive its dispersion relation. With some differ¬ 
ences, the spectrum is similar to that found for EM wave propagation in covalent crystals. 
We then add to this continuum theory isolated defects, so as to obtain a more realistic 
model of defects embedded in a frequency dependent dielectric medium. We do this both 
with a single isolated defect and with an array of isolated defects, and so compute how the 
S-matrices and the band structure change in a dynamic medium. 
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1. Introduction 


In this paper we examine photons interacting with harmonic oscillator impnrities in an 
optical hber. Becanse this problem is inherently two dimensional, we are able to attack it by 
applying or extending varions techniqnes from two dimensional defect and bonndary held 
theory ([Q 0 0). This permits ns to compnte exactly qnantities previonsly inaccessible. 
Here, for example, using methods akin to thermodynamic Bethe ansatz (TBA), we compute 
thermodynamic potentials. And through extending the notion of boundary states, we 
exactly compute correlators. 

Normally the impurity atoms in optical problems are modeled as two-level systems. 
Although this makes the problem non-linear, some success has been had in hnding exact 
solutions in two dimensions. Under the rotating wave approximation, a continuum array 
of defects has been shown to be integrable (0). There, using algebraic Bethe ansatz, the 
spectrum and the S-matrices of the theory were found. In it was shown that a single 
two-level atom can be mapped onto a Kondo impurity. With this mapping they were able 
to compute thermodynamic potentials exactly and to obtain form factor expansions of 
some correlators. But such successes are limited. Because of the theory’s non-linearity, 
integrability is key to obtaining exact solutions. However integrability limits the type of 
defect conhgurations one can consider. Outside the rotating wave approximation, an array 
of two-level atoms is not known to be integrable. To then say something one must resort 
to perturbation theory (see |P]). And although a single two-level atom can be treated, 
two or some larger number cannot. Nor does integrability guarantee that all interesting 
physical quantities can be accessed. Photon correlators for two-level systems are still 
only known approximately. And the computation of correlators involving the impurity 
degrees of freedom remains an open problem. To sidestep these difficulties we choose to 
model the atoms as harmonic oscillators. In doing this we cease to be limited by particular 
conhgurations of atoms under certain approximations and we gain access to a wider variety 
of quantities, including both exact photon and impurity correlators. 

The replacement of two-level atoms by harmonic oscillators is not necessarily a bad 
one. Providing the intensity of the photon held in the hber is low, only the oscillator’s 
ground and hrst excited state are likely to be occupied. Under these conditions, the 
oscillator mimics a two-level system. We will see later, for example, how the spontaneous 
emission spectrum mimics that of a two-level atom. However, we cannot hope to model 
such phenomena as strong held resonance huorescence. 
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We exploit our ability to study arbitrary configurations of defects by examining both 
periodic and continuum arrays of defects. With a periodic array of defects, a band structure 
of photons naturally appears. This band structure is equipped with gaps of two types: one 
arising from being on resonance and a more typical one of gaps appearing at the edges of 
the Brillouin zones. 

In the case of the continuum array of defects, we derive a dispersion relation for the 
photons propagating among the defects. We find that the dressed photons are described by 
the typical polariton dispersion relation. Again exploiting our configurational fiexibility, 
we examine theories where isolated defects are added to the continuum array. We look at 
two examples: an isolated defect and a periodic array of defects. With the isolated defect 
we compute the scattering matrices of the polaritons oh the defect. In the case of the 
periodic array, we calculate the band structure obeyed by the polaritons. 

The ability to consider defects living in a dynamic medium lays the ground for the 
consideration of photon localization. By choosing the energy of an isolated defect to lie 
within the gap between the two polariton branches, there is some chance to find a photon 
in a bound state. We will consider this possibility in an upcoming work. 

The paper is organized as follows. In the first section an overview of the model as well 
as its reduction to two dimensions is given. In the next section the scattering matrices are 
derived, followed by a section giving the TBA analysis of the thermodynamic potentials. 
In the next section correlators are computed. There we compute both time ordered and 
non-time ordered correlators. The non-time ordered correlators, important as they often 
represent what can be physically measured, are accessed via analytic continuations of 
temperature correlators. In the next section we consider multiple defects and in particular 
periodic arrays of defects. In the final two sections, we consider a continuum array of 
defects, first alone, then together with embedded impurities. 


2. Overview of Model 

In this section we describe the coupling between a single harmonic oscillator defect 
(placed at the origin) and a photon field in the environs of an optical fiber. The oscillations 
of the defect have two possible physical sources. The first, going under the name “atomic 
polarizability”, arises from the distortion of charge within an atom. One can crudely model 
this by representing the atom as a charged shell tied to an immobile nucleus through a 
spring of strength K = mujQ (m is the mass of the shell and equal to Znie where Z is 
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the atomic number and mg is the mass of an electron). The shell then executes simple 
harmonic motion around the hxed nucleus. Typically the energy scale of these oscillations 
is huo lOeV. 

The second source comes with the name, “displacement polarizability”. Here the 
defect is considered as a pair of oppositely charged ions that oscillate about their centre 
of mass. Because these oscillations are phonon-like (i.e. the heavier ions are oscillating 
and not the electronic shell), their energy scale Hloq is typically 10“^ ~ 10“^eV. In this 
section we will consider the defect to be of this latter type. At energy scales of the former, 
we would also need to consider the atomic polarizability of the hber medium in which the 
photons propagate and in which the defects are embedded. Instead we are able to treat 
the medium of propagation as vacuum (or equally easily, a medium of static dielectric 
constant, cq). We will delay the inclusion of the dynamic properties of the medium to 
Section 7. 

The oscillator, minimally coupled to the photon held in the radiation gauge, has the 
Hamiltonian 


77= — + 

2m 


mu, 


0q2 


^A(x = 0) + ^ 

2m 2mc^ 


A2(x = 0). 


( 2 . 1 ) 


By now introducing, in the standard way, the variables 


^ -:T7yP; 

^ / {2mhuo)^ 

1=[S-,S+], 


this Hamiltonian can be recast as 

(S+S-+ 1/2) - - (^) (S++S-)A(x = 0) + ^A"(i = 0). (2,3) 

The interaction terms are then of two sorts: one coupling the photon held to the oscillator 
and one giving the photon a “mass” at a: = 0. Throughout our analysis we will keep this 
latter term. The term contributes equally to the quantities which we will be interested 
in calculating. Its omission, on the other hand, thereby violating gauge invariance, leads 
to unphysical results. Fortunately, its inclusion does not overly complicate the analysis. 

As we ultimately want to derive equations of motion, we want to write the above as an 
action. Because S~^ and S~ can be thought of as conjugate variables, i.e. [ihS~, S'"*"] = ih, 
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this can be easily done. The action is 


Ssho = 


j dt {-ihS-dtS^ - H) 

j dt {-ins-dtS+ - huoS+s 

1/2 


(2.4) 


+ 


(5" + S-) A(. ^ 0) 


2mc^ 


A^{x = 0 ) 


In writing this we have dropped the zero point energy of the oscillator. 

We now focns onr attention on the free dynamics of the photons. The action deter¬ 
mining these dynamics is 

= (2.5) 

For a hber, Sa can be dimensionally reduced. We begin by setting the scalar potential Aq 
to zero. We then take the hber along the T—direction. A is necessarily in some transverse 
direction. Assuming only one direction of A is excited, say h where h-x = 0, we can write 
A = An. If A is independent of the transverse directions, y and z, Sa reduces to 

= ^ J dtdx {c~‘^{dtAf - {dj:Af) , (2.6) 

where D is the cross-sectional area of the hber arising from performing the integral J dydz — 
D. We can rescale A to turn Sa into an action describing a 2-d free massless boson: 


A 


Di/2’ 


(2.7) 


SO that 

= j dtdx . (2.8) 

We list factors of c and h so as to be able to easily provide estimates of various couplings. 
The full action we consider is then 


S 


Sa + S 


sho 


j dtdx [c ^(9t(/))^ — (9x</>)^) + J dt ^ — ihS dtS^ 
-hujoS^s- + {nc)^/^x {s+ + s-)(f)- j. 


where A and are 


/ ^0 e_ 

V 2mcD ) c ’ 

e^h 

2mc^D 


(2.9) 


( 2 . 10 ) 
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has units of mass and A has units of inverse length. Because the bound, approximately 
by 10“^, on the hber diameter, D, is large compared to atomic length scales, both A and 
are small compared to the scale set by hug ~ 10“^ev, the typical oscillator energy. We 
have 


Xhc 

hUg 

A^C^ 

hUg 


10 


-7 


( 2 . 11 ) 


10 


-14 


We thus see that the mass perturbation A^ is on the same order as A^, as required by 
gauge invariance. When h = c = 1 the two scales are governed by 


A^ = A'^UJg. 


( 2 . 12 ) 


Thus A^ cannot be ignored in favour of A, the photon-oscillator coupling. Because of 
(2.12), we will express all derived quantities such as scattering matrices, band structure, 
correlators, etc., in terms of A^. 

Though both couplings are small compared to hug, the perturbations are still able to 
lead to interesting physics. It is possible to introduce other length scales into the problem 
by considering multiple defects, which in combination with the existing scales can lead to 
observable effects. We will consider such cases in Sections 6, 7, and 8. But for now we 
continue to develop the single defect formalism. 


3. Calculation of Defect S-Matrices 

Generally given a theory with particles A{k), an integrable defect perturbation pro¬ 
duces scattering via the defect S-matrices, T{k) and R{k), as follows: 

BA+ik) = T{k)A-{k)'D + R{k)BA+{-k), (3.1) 


where D is the defect operator. D is analogous to the boundary operator B discussed in 
0. A± represents particles who live on the right {x > 0)/left (x < 0) side of the defect. 
T{k) is then just the amplitude that the particle, A_|_(/c), makes it across the defect, and 
R{k), the amplitude that it reflects back. 

To derive T and R we reconsider the action in (2.9). Setting h — c — 1 and trans¬ 
forming to Euclidean time, t —it, this becomes 


Stt 


S =:^ I dtdx {{dt(j))‘^ + {da:(f))‘^) 

+ [ dt ujgS+S- - S-dtS+ - X{S- + S+)(j) + 


(3.2) 
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By varying S , S'"*", and (/>, we obtain the following eqnations of motion: 

0 = Xd{x){S~ + S~^) + — — 2d{x)A'^(f)-, 

0 = 5(a:) {-X(P + LOoS^ T dtS^) . 

The second of these two has the steady-state solntion 

:± _ \ /'°° 


= X 


dke 


iVoTk' 


(3.3) 


(3.4) 


where (f){k) is dehned throngh 


(j){x = 0) = / dke^^'4>{k). 


(3.5) 


We can thns eliminate from the hrst eqnation of motion to obtain 

0 = 2ljoX‘^S{x) [ dke^^ ^ + dx'f) ~ 2A^(j)5{x) 


To solve this eqnation we write cf) as 


(3.6) 


4>{x) — 4>j^{x)6{x) -|- 4>-{x)6{—x)^ 


(3.7) 


where (^± are free massless bosons with mode expansions 

1 


^±{x) = 




yl±(fc)e-*l^l+^^^ + yl^(A:)e‘l^l 


—ixk 


(3.8) 


We take (l){x = 0) to be dehned as (/)(0) = l/2((/)_|_(0) -t- (/)_(0)). So onr eqnation of motion 
can be written as 


0 — cc’O'^ 


dke 


kt4>+ I 7 /„ / O , A 2/ 


lVq — k‘^ 47r 


+ — - dx<l>-m - A"((/>+(0) + (/>_(0)). (3.9) 


Second order derivatives have vanished becanse {d^ -|- d‘^)(f)± = 0. 

To derive T{k) and R{k), we divide the above eqnation into its ‘-t-’ and parts, 
say /(A_|_(/c)) -|- g{A-{k)) = 0, and so interpret it as vanishing when acting on the defect 
operator in the following manner: 


JDfiA+ik))+g{A.ik))JD = 0. 


(3.10) 
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Note carefully the orderings of D and the ’s and ’s. Substituting the mode expansions 
into this equation together with the continuity equation 

D(/>+(0) = (/)_(0)D, (3.11) 


hxes T{k) and R{k) to be 


T{k) 

R{k) 


ijk^ -‘^o) 

i{k‘^ — Uq) — AnA'^k ’ 
AnA'^k 

i{k'^ — uJq) — AnA^k' 


(3.12) 


As a check on the validity of this result, we see T{k) and R{k) satisfy the defect unitarity 
conditions: 

T{k)R{-k) + T{-k)R{k) = 0; 

(3.13) 

R{k)R{-k)+T{k)T{-k) = 1. 

These conditions are easily derived through applying ( p.l| ) twice, i.e. allowing the particles 
to scatter through the defect twice. 

The analytic structure of T{k) and R{k) is unsurprising. Both have poles at 

kp = ±wo(l - f )^/^ - i27rA2 (3.14) 

Uq 

The real part of kp describes the energy splitting between two oscillator levels. This 
splitting has been shifted slightly from ujq in a process analogous to the Lamb shift. Minus 
the imaginary part of kp equals the inverse lifetime of an excited state. We then have 


T ^ — —Imkp = 27rA^. 


(3.15) 


This lifetime is what is similar to what is found in perturbative treatments using the 
master equation (see for example [0). This is unsurprising. If we were to subject the 
theory to a perturbative analysis we would hnd that the photon self-energy (i.e. the 
vacuum polarization) consists of a single diagram of lowest order and so should have 
no contributions higher than A^. A similar situation is found in perturbative analyses of 
the simple defect theories in [Q. 

The interpretation of this pole differs from the standard one in integrable Euclidean 
QFT. Normally the imaginary part of kp would be interpreted as related to the energy of 
an excited oscillator state and the real part as the decay time of the excited state. But 
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here the roles have been reversed, as if we still were in Minkowski space. This standard 
interpretation arises from assnming that the pole in the S-matrix arises from an effective 
propagator of the form Bnt here this is not the case. In Minkowski space, 

the pole appears as ~ (/c^ — Uq)~^. Both k and uq are energies. In going to Enclidean 
space, we then have k —ik and uq —iojq. Thns the pole structure retains its form in 
Euclidean space. 

4. Thermodynamics 

Free photons in a cavity at a temperature T are governed by the well-known Planck 
distribution: 

dk 1 

= ^4.1) 

The mean energy per unit length is given by 

/ C)0 

dk\k\pik) = -T\ (4.2) 

-oo 6 

The above expression for U{T) is the one-dimensional analog of the Stephan-Boltzman 
law. U (T) can also be expressed as 

U{T) = ^dT\ogZ, (4.3) 

where L is the length of the cavity, and 

T / \ 

logZ =-/ dA;log(l-e-l^l/^), (4.4) 

27r V / 

is the free energy of the non-interacting photons. We will now compute the modihcations 
of the above blackbody theory due to the presence of the impurity. We place the impurity 
at the origin at x = 0, let the cavity have length L, —L/2 < x < Ll2, and further impose 
periodic boundary conditions at the ends of the cavity. 

For a single impurity, we can fold the system so that it is formulated as a theory on 
the half line with the impurity on the boundary at x = 0. This kind of folding was used 
in the work Q in their study of edge states in the z/ = 1/3 quantum Hall system. Let 
4> — 4>l + (pRi where +1) and (/»/j(x — t) are the left and right-moving components of 
the massless scalar held in the bulk. Dehne even and odd helds as follows: 

t) = t) + (j)R{-x, t)) / a / 2 ; 

= {(i)L{x,t) - (j)R{-X,t)) /^/2■, 

(^%{x,t) = (j)l{-x,ty, (i)R{x,t) = -(j)l{-x,t). 

8 


(4.5) 



Also define = 
half-line x < 0: 




+ 


^R- 


Then the action (^) can be rewritten as a theory on the 


S = 


Stt 


dt / dx 


' —oo ^ —oo 


/ oo / \ A 2 \ 

dti- iS-dtS+- u;oS+S- + ^{S+ + S-) r - W')d • 


(4.6) 


Note that (j)^ and (j)° deconple, and only conples to the impnrity. 

The interaction of the even photons with the impurity at the boundary are now 
summarized by a reflection matrix TZf,- The algebra (|3.1|) now only has a reflection piece: 


BA{k) = TZe{k)BA{-k), 


(4.7) 


where B is a boundary operator. A similar computation to the one in section 3 gives 

i{k'^ 


n,{k) = 


Uq) + AttA^/c 
i{kP‘ — lOq) — AttA^/c 


(4.8) 


A similar reflection matrix describes the odd photons. However as the odd photons do not 
couple to the defect, it is trivial, equaling unity. 

Thermodynamic properties can all be expressed in terms of TZe- The quantization 
condition on the moment um k of an even photon on the half-line of length L/2 is 


ne(k) = 1, k>0. 


(4.9) 


We need only consider k > 0 since it is implicit that the wave-function satisfying the above 
quantization condition contains both a right-moving and reflected left-moving particle. 
Taking the /c-derivative of the log of ( |4.9|) yields 


1 


Pe{k) = — 


271 27rL 


5tlogKe(fc), 


(4.10) 


where Lpf,{k)dk is the number of allowed states between k and k + dk. Introducing a 
density of occupied levels Pe{k), one has for the partition function of the even held 


logZe = + Se, 


(4.11) 


where Ue = L dkkpe{k) and Se is the entropy: 


Se=L 


dk{pe -F Pe) log(pe + Pe) “ Pe log pe “ Pe log pe 


(4.12) 
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Minimizing logZ^ with respect to pe, one finds that the Planck distribution (ITT) for the 
even photons is modified to 


Pe{k)dk = t/T_i (4 - \ dk, 


(4.13) 


and the even free energy becomes 

log Ze = -L 


Pe(k) log (l - e • 


(4.14) 


We now put this all together to find a total energy for the system. 

The mean energy density separates into bulk and impurity contributions: 

[/(r) = [/buik(T) + [/i^p(r). 

b^buik is a sum of equal contributions from the even and odd sectors of the theory, 

k 1 

t^bulk = + t^bulk = 2 X / dk 

Jo 


(4.15) 


271 e^!'^ — 1 


(4.16) 


Hence tfbuik = 17imp is the piece from the even sector that scales as L ^. Hence 


4A^ k 

t^imp(T) = ^ I dk 


k'^ + UJn 


L Jo — 1 (/c^ — Uq^ + (47rA2/c)2 


(4.17) 


In the limit A —> 0, the photons are decoupled from the impurity, and one expects LUiuip 
to be just the energy of a decoupled harmonic oscillator: 


hm = T^dr log ( j;] 


1 


— ^0 


n=0 


qUJoIT — 


(4.18) 


We have set the ground state energy of the harmonic oscillator to zero since in the scattering 


theory description the vacuum is defined to have zero energy. One can confirm that (|4.17| ) 
has the proper zero coupling limit, using the identity: 


hm 


1 


1 


Then, 


aZ/ujo^o {kP- — u>q)^ + (47rA2/c)2 SA'^Uq 

|im Lf/i„p(r) = ojo— 

which agrees with ([4.18|). 


d{k - (4;o)- 


^ujQ IT — 2 ^ 


(4.19) 


(4.20) 
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In common physical situations, Uq S> (47rA^)^. Making this approximation in ([4.17|) , 
one can extract the following low and high temperature behaviours: 

27r^ 

— — , T<wo; 

6 1^ Uq ^^ 21 ) 

UimpiT) T > Uq. 

5. Computation of Correlators 

In this section we will consider how to calculate several types of correlators that are 
useful in quantum optics problems. We begin with time-ordered correlators. 


5.1. Time-Ordered Correlators 


There are two possible pictures in which to view the defect. The hrst, as in section 3, 
treats the defect as a point in space. Particles then scatter to and from the defect. Because 
the Hamiltonian of the theory is altered by the defect on any constant time surface, the 
particle spectrum is altered. This is reflected in the linear relations between the particle 
creation/destraction operators as in 3.1 and 3.2. In the second picture the defect is a point 
in time, an initial condition. Away from the defect the Hamiltonian of the theory remains 
unchanged and so the particle spectrum in the theory is the same. This is of particular 
advantage in calculating correlators. The defect, as an operator, is expressible in terms of 
the unperturbed Hilbert space and the correlator is then calculated by inserting the defect 
within the free helds forming the correlator in a time-ordered fashion: 


{T{Oi{xi) ■ ■■On{Xn))) 


(D) 


(5.1) 


On the l.h.s. of the above equation, the helds are the full Heisenberg helds whereas the 
time evolution of the helds on the r.h.s. is governed by the free Hamiltonian. D (to be 
distinguished from the D of section 3) then is no more than the S-matrix for the theory. 
It has a particularly simple form 


D = exp 


dkT(k)AUk)A_(k) + / (-i) + >l_(i).4_(-i)) 


(5.2) 


The ± subscripts on the photon operators indicate on which side of the defect the particle 
lives: -|- for t > 0 and — for t < 0. The hrst term in the exponential transmits parti¬ 
cles through the defect, destroying a — particle and replacing it with its -f counterpart. 
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Likewise the latter two terms represent reflections off the defect. Given D’s form, it is com¬ 
pletely represented in terms of one particle scattering amplitudes. This is a consequence 
of the theory’s integrability. The form of the terms in the exponential are constrained by 
momentum conservation; the total momentum must sum to zero. Because all the terms in 
the exponential commute with one another, there are no ordering problems in representing 
D. We have assumed here that (0|D|0) = 1, that is the defect entropy is zero. With finite 
temperature this will not the case. 

To distinguish between the two pictures, space-time coordinates in the first picture 
(the physical picture with the defect as a boundary condition) will be written as (x, t) and 
the space-time coordinates in the second picture (with the defect as an initial condition) 
by (x, t). The two are then related by (t, x) — {x't). To compute the elements of D, T{k) 
and R{k), one first expresses the action in the rotated picture, obtained from 3.3 via a 
Euclidean rotation and a reflection: x ^ t and t —x. As before we can then derive an 
equation of motion with the eliminated: 

0 = ^ (Sri>+( 0 ) - ar<^-( 0 )) - a"(,^+(o) + ^_(o)), (5.3) 

where (j)+ and (f)- are defined via cf) — -|- As in section 3, we have a 

continuity equation, 

=0) ^ =0). (5.4) 

To interpret the action of the above two equations on D, we divide each into a part 
involving (j)^ and a part involving (j)-. So both have the form /_(_((/)_(_) — /_((/>_) = 0. In 
terms of D these become 


/+(</'+)D^(</'-)|0) = (0|/i((/>+)D/_ ((/>_). (5.5) 


g{4>-) and h{4>^) represent arbitrary insertions of the fields cf)- and between D and 
the vacuum state. (|5.5|) tells us that as the fields /_|_((/)_|_) cross the defect D, they are 
transformed into /_((/>_). 

To determine the form of T{k) and R{k), we expand out D to first order in T{k) and 
R{k) and substitute into the above two equations. Consistency at the one particle level 
then yields 


Tik) 

R{k) 


Uq -H /c^ 

u>Q + k‘^+ 4:7rA^\k\ ’ 
— 47rA^|/c| 
iVQ + k'^+ 47rA2|/c| 


(5.6) 
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T{k) and R{k) can be related to the scattering matrices of section 3 via an analytic 
continnation. We see 

f{k){\k\^-ik,k^i\k\)=T{k), (5.7) 

and similarly with R{k). This analytic continnation is effectively a rotation of the defect 
back to the t-axis. Here energies and momenta are continned via E —> —ip and p iE. 

Knowing D we can now go on to calcnlate correlators. We will compnte the two-point 
fnnction of the photon held. We have for this correlator, 


(T((/)(T,f)(/)(T',f))) =e(t - t) 


9{t)0{t ){(l){x,t)(l){x',t )D) 

+ 9{t)d{—t ){4>{x,t)'D(l){x',t )) 
+ 9{t)0{—t ){D(l){x,t)(l){x',t )) 
+ [t ,x^x) 


(5.8) 


Each of the pieces of {T{(f)(f))) is easily calcnlated as (/>(/) only conples to the two particle 
contribntion to D. So 


{(j){x, t)(j){x'X)'D) 
{(j){x, t)U(j){x'X)) 


{4>{x,t)4>{x',t ))free + 2 f dk ^[^^ 6 *■*"*“* cosikCx — x')); 

k 


\k\ 


{D(f){x,t)(f){x')) — {(f){x,t)(f){x' ,t'))fT:ee +‘^ [ dk ^^^^ 6^^ COs(/c(T — T')). 

Jo ^ 

(5.9) 

The correlator {T{(^4>)) has now been calculated in the picture with the defect as initial 
condition. To obtain this same correlator in the physical picture, we rotate and rehect the 
system back via T —> t, t x. If, for example, we supposed a: = x' > 0 in the original 
system, the correlator becomes 


{T{(i){xR)(i){xR')) = ((/)(T, t )(/)(T', t )D) 

= x)(f){t', x))free + 2 / cos{k{t - t')). 

Jo ^ 


(5.10) 


The correlator only depends on R{k) as we have assumed that both space points are on 
the same side of the defect. 
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5.2. Non-Time Ordered Correlators 


So far we have only considered time-ordered correlators. Bnt correlators not so ordered 
can also be easily compnted, a good thing as the correlators in which one is often interested 
for qnantnm optic calculations are not time-ordered. For example, it common to want to 
compute Glauber correlators, correlators of normal ordered products of the creation and 
destruction pieces of the photon helds. In this section we will consider three examples 
of non-time ordered correlators: the photon density in the hber at hnite temperature, 
the spontaneous emission spectrum for an excited oscillator state, and the correlators 
of the defect degrees of freedom. We will do these calculations in a hnite temperature 
formalism. It turns out the calculations are easier to do in such a framework. Euclidean 
correlators which are not time ordered are not convergent for all of spacetime. This makes 
the rotation of the defect between viewing it as a boundary condition and viewing it as 
an initial condition ill-dehned for such correlators. But in a hnite temperature formalism, 
unordered correlators can be related to time-ordered temperature correlators for which the 
rotation, {x,t) ^ poses no difficulties. 

Temperature correlators are computed in much the same ma nn er as the time-ordered 
correlators of the previous section. The primary diherence lies in the mode expansions. 
As we are at hnite temperature, T, we expand the helds on a cylinder of circumference 
(3^1/T = 47 r/: 

= (t + ix)/2; 

n>l 

n>l 


We have divided (f) into its left and right moving components. Here the mode expansions 
are periodic in x (and not t) as we intend to eventually interchange time and space (so as 
to return to the original system). 

The defect operator D for T 7 ^ 0 is constructed in exactly the same fashion as in the 
T = 0 case. We substitute the mode expansions into the equations of motion and hud 


0 — —{dt(t)Jr ~ dt(p-) — A^((/)+ + (p-) + X{S^ + S ); 




i>i 


n 


1/2 


Cc^Q —p 


-{4>-n + 4>n) + 




-{(pn + (p-n) 


(5.12) 
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where the Matsubara frequencies are given by cOn = n/2l. From this we construct D as 
before: 


T) = g exp 


T{Un){(l)tr,(l)n + <P-n<P-u) + + </>n </'n ) , 

n>l 


(5.13) 


where T and R are as in (5.6). In the T = 0 case we were able to set g, the boundary 
entropy to 1. Here we cannot do so. The partition function of the system is 


Z = (0|e‘‘''De‘»''|0), 


(5.14) 


where time extends to some 1 1 and 1 2 on either side of the defect, {t i + 1 2 thus dehnes 
the length of the cylinder which we are on.) As 1 1 , t 2 —^ 0, we expect Z to reduce to its 
defect contribution. Thus 


Zimp = _ lim Z = (0|D|0) = g. 

t i,t 2 —>0 


(5.15) 


Our TEA calculation from the previous section gives us 


logs' 


-4A^ 



(uj^ + P)\og{l-e ^/'^) 
{uq — / c ^)^ + 


(5.16) 


We point out that knowledge of g is not necessary for the computation of correlators as it 
is a constant which is always normalized away. 

We now go on to compute the number of photons, n{uj), in the hber. n{uj) is given in 
terms of the two point photon correlator: 


n{uj) = — / dte (5-17) 

Stt J 

where are the creation/destruction pieces of the held. These pieces are dehned for 
non-free held theories to be (see [^) 

= ±^. j dt'ZppR. (5.18) 

The factor c<;/27r in (5.17) serves to remove the relativistic normalization present in the 
mode expansion of the photon helds. is not time ordered and as indicated previously 
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would be difficult to compute directly. However by the fluctuation-dissipation theorem we 
can reexpress this correlator in terms of a retarded Green’s function: 


dte {(jf{x^t)4>{x^Q)) = 


ol3uj _ 


lmGR{-u)6{u!)- 


(5.19) 




This sort of identity between correlators is easily proven by an explicit evaluation of the 
thermal trace together with an insertion of states between the operators (see 0). The 
advantage of writing n(c<;) in this form is that Gr{uj) is available directly as the analytic 
continuation of a temperature Green’s function, 


Gr{uj) — Gri—ii^ + ^) — 


e'-‘^^^(j){x,t)(j){x,0))': 


(5.20) 




and Gt is exactly what we are set up to calculate. 

To compute Gr = —{(f){xG)4>{^^ 0))t we go into the rotated system. Assuming x > 0 

we have ^ _ 

Gt = -(0(x,t)(/)(x,O))i = -((/)+(x, ?)(/)+(x,0)D) 

= -7 — cos(wnt)[l + (^-21) 


n>l 




Following the prescription above, we then And for n{Lv), 


n{u) 


2e{uj) 

- 1 


1 -|- sin(2ci;x) 
— cos(2ci;x) 


47rA^Ci;(ct;^ — Uq) 

(w^ — Wg)^ + IAtt^A^w^ 


(5.22) 


We see that we obtain the Planck distribution (times two - one for right moving photons 
and one for left moving photons) plus a correction term. This distribution should not 
be compared with the TEA calculation of the previous section. The distribution there is 
not that of pure photons but rather of their even and odd combinations that have been 
averaged over the length of the system. 

We now move on to calculating the spontaneous emission spectrum, P{uj), arising 
from the decay of the oscillator from its first excited state. As the electric held is given in 
terms of (/» by i? cx dtcj)^ we define the power spectrum to be 


P(ul) 


1 


4jj-2 


dtdt' ^(5'“(0)(/)(x, t))((/)(x,t')5'"'"(0)). 


(5.23) 


lA 












In the second line the fonr-point correlator has hrst been factored into two two-point 
correlators and replaced with (p, both possible as only creates or destroys single 
particles, and then the resnlting expression has been integrated by parts. These steps are 
only possible at T = 0. However we will let T 7 ^ 0 in order to compnte the two-point 
fnnctions and then take T —> 0 at the end. We can make a farther simplihcation by noting 


that 


J die ( 0 )(/)(a:, t)) 




* 


(6.24) 


We thns need only to compnte {(p{x,t)S~^{0)). 

As before we compnte {(p{x,t)S~^{0)) throngh a temperatnre correlator. In this case 
we have 


J ^ ^ -e-i 3 u {Gri-iuJ + e) - - e))- 

r(3 

GtM = - dtG^-^{<P{x,t)S+{0))T. 

Jo 


(5.25) 


Onr expression involves both the retarded Gt{—x> -t- e) and the advanced Gy (—a; — e) 
Green’s fnnctions as the combination (pS^ is not Hermitean. Again this identity can be 
verihed by evalnating the thermal trace and inserting a complete set of states between p 
and S~^. The temperatnre correlator is compnted as before by rotating the system. Using 
(5.12) we then hnd, 


{p{x,t)S+{0))T = ((/)+(T,f)*S+(0)D) 


-E 

21 ^ 


T{Ur 


n>l 


Ur 


Uq — iUn Uq -|- iUn 


(5.26) 


The Matsnbara decomposition of this correlator is 

GtM = - e{-Urr)e--^f{-Urr)). (5.27) 

Pntting everying thing together and taking T —0, we hnd the power spectrnm to be 
2A'^uo6(u) {uq — w^)^ 


P{u) = 


{uo - up 


D 


1 + 


COs(2wa:),, O , . 2 n 2 -,^_2a4, . 2 ^ 


D 

sin( 2 c<;a:) 

D 


-{{uq—u^) — 167r^A"^c<;^) 


(87rA^ct;(c<;Q — u^)) 


(5.28) 


D = (uq — u'^Y + IBtt^A" 


u . 
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Typically the power spectra is computed for near resonance (w ~ wq)- Near resonance 
P{uj) reduces to 

4A^uJo sin^(ci;a:) 


P{u) = 6{uj)- 


(5.29) 


(wo — 

We see that P{u:) takes on the familiar Lorentzian lineshape with halfwidth equal to 47rA^. 
Near resonance the lowest order term of P{uj) has the form predicted by classical theory 


(see [n|) and the same form predicted by hrst order perturbation theory employing the 
master equation (see [0). The spatial dependence in the power spectra is an artifact of 
not including the power from the magnetic held. Doing so leaves P{u:) independent of x. 

This basic Lorentzian lineshape is also shared by a two-level system interacting with 
a photon held near resonance (see 0]). Where the physics of the the two-level system 
disagrees with that of an SHO defect is in the presence of an external driving held. For a 
SHO defect the response spectrum resulting from an external held is a h-function. This is 
only true for a two-level system for a weak driving held. As the driving held intensity is 
increased, the 5-function spectrum of the two-level system broadens and then splits into 
three Lorentzians, one centered at the driving frequency and two sidebands equidistant 
from the central one. The two systems diverge in behaviour in large external helds as the 
two-level system saturates in its upper level for such held whereas the SHO can keep being 
excited into ever higher levels. At low driving intensities the SHO spends its time mostly 
its in ground state or its hrst excited level and so mimics a two-level system. 

As the hnal illustration of computing non-time ordered correlators, we compute the 
correlators of the defect degrees of freedom, i.e. (S">'(t)S''^(0)), 7,5 = ±. Even though the 
defect degrees of freedom have been integrated out, we can still compute their correlators 
as we know how are related to the photon held. As before {S'^(t)(0)) is related to a 
temperature correlator via 




dt 




9{u! 


2711 — e 
rP 

iuJn t i 


+ e) - - e)); 


GtM = - 


{S'^{t)S^{0)) is computed by rotating the system 

{S^{t)S\0))T ^ {Sl{x)BSt{0)) 


T- 


(5.30) 


(5.31) 


As and are dehned at x = 0, we have a choice as to their ordering with respect to 
D. The expressions for the correlators are independent of this choice (as they should be). 
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We have chosen in the above to place and on either side of D. Using (5.12), the 
correlator on the r.h.s. of the above is easily calcnlated: 


(^7(x)D^i(0))T = 


^icJn t 


^ — t 


+ 


i>] 


n [(wo - + iSuJn) + i'yUn ){(^0 - iSuJn), 

(5.32) 

Then by ( b.30|) , {S^S^){uj) eqnals 


{S^S^){u) = 


2A^uoO(uj) {uJq — w^)" 


u{l — e ^^) (wo + 7‘^)(‘^o — 5u) {uq — 


(5.33) 


One shonld notice that the correlators {S~^S~), {S~^S~^), and {S~S~) do not vanish. This 
is a resnlt of the conpling between the photon held and oscillator. Becanse the interaction 
is {S~^ + *S'“)(/) (and not + S ~and S~ have both creation and destrnction 

parts and so undergo vacuum huctuations. 


6. Multiple Defects 

The analysis of a single defect in the previous sections can be easily extended to 
multiple defects. We will hrst consider a set of n defects placed arbitrarily and then move 
on to a periodic array of defects. 


6.1. n Defects 

The action for a set of n defects is 


S 


dtdx 


1- {(Mf + ipAf) 


+ 5^ - Xi) (-sra,s+ + uJiSfS- 

i=l 


Xi(S- + S+)4> + 


( 6 . 1 ) 


where the i-th defect has been placed at XiJ^i > Xi> Xi-i. We have allowed each defect 
to have different coupling strengths to the photon held. The equations of motion for this 
action are 


1 

0 = ^h(x - Xi) [XiiSr +S+)- 2A2(/)] + — {df(f + dl(i)) ; 

i—1 

0 = 5{x-Xi) + uOiSf" =F dtSf). 
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By assuming cj) has a free solution in the intervals in between the defects, 


4) = - Xn) + - x)6{x - X^-l) H- 

+ (j)2d{x2 — x)9{x — xi) + ~ ^)-i (6-3) 

0 = + dl)(t)i, 1 < i < n + 1, 


the above equations of motion can be reduced to 


It - n 

0 = - Xi) uJiX^ i 

i=i L 


dke^ 


where (j)i is, as before, 


- w? 


+ + 4'i) ~ [dx^'i+l — dx^i) 

47r 


I dke^^4}{k,Xi). 


(6.4) 


(6.5) 


We have thus obtained a set of n uncoupled equations. These can be solved in same fashion 
as in Section 3. 

Let Di be the defect operator for the i-th defect. We dehne the transmission and 
reflection matrices across this defect by 


D..4.+i(fc) ^ ri.4i(fc)Di + fii(fc)Di.4i+i(-fc), 


( 6 . 6 ) 


where the Ai(/c)’s arise from the mode expansion of (pi- ( |6-4|) then hxes Tj and Ri to be 


Ti{k) = T{k,Ui, Ai), 


R^{k) 


— 2ikx 


^R{k,Ui,Ai), 


(6.7) 


where T(/c, uji, Ai) and i?(/c, Ui, Ai) are the matrices in (|3.12|) . We see that Ri has picked 
up a phase, relative to R, a result of a reflection away from x = 0. 

Ti and Ri describe the transmission and reflection across the defect at Xi. However 
the transmission and reflection amplitudes across all n defects is considerably more com¬ 
plicated; one has to consider multiple reflections and transmissions as the photon traverses 
the set of defects. To simplify the computation, we introduce a transmission matrix dehned 
via 

^ / Ai{k) 

A,+^{-k)) ~ * VA*(-^) 


( 6 . 8 ) 
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Using (|6.6|) and its conjugate, 


Ai(-fc)Di = Ti{k)'DiAi+i{-k) + Ri{k)Ai{k)'Di, 


M, is found to be 


T-\-k) R^ik)Tr\k) 

R,{-k)T-\-k) T-\k) 


We can then write down the transmission matrix across all the defects: 


(6.9) 


( 6 . 10 ) 


Mt = ■ Ml (6.11) 

The overall scattering amplitudes, Tt and Rt are dehned by 


R^TAn+iik) = TtAi(/c)Dt + RT{k)Y)TAn+i{—k), (6.12) 


where Dt = ■ ■ • Di. These amplitudes are then related to Mt via 


IVI'p — 


T^\-k) 

RT{-k)T^^{-k) 


RT{k)T^\k) 

T-\k) 


(6.13) 


As it stands Mt is not computable analytically (though certainly numerically) for an 
arbitrary set of n x^’s, A^’s, and w^’s. As such we will only consider a particularly simple 
case: the case of two defects. Taking xi = 0, TTik) and Rrik) can the be computed to be 

= 1 L = Ti(k)T2(k) f^(R2{k)R2{k)r 

l-Ri{k)R2{k) ^ 

oo (6-14) 

RT(k) = R 2 (k) + , = R 2 {k) + Ri(k)n(k) f^(R 2 (k)R 2 (k)r. 

l-Ri[k)R2{k) ^ 

The latter equality in the above two equations shows that Tt and Rt can be both computed 
through adding up multiple reflections and transmissions as if this problem was a two 
interface problem in classical optics. 

As with T(/c) and R{k) in Section 3, both Tt and Rt exhibit poles as a result of the 
oscillator defects. Rt has simple poles at 

kp = ±u;2{l - ^)^/2-z27rA2, (6.15) 

^2 
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indicative of reflection off the second defect. There is no similar pole for the first defect as 
we are scattering from the right. If we were scattering from the left, a pole in the variables 
oji and Ai wonld appear in Rt- Rt and Tt have another set of poles dependent npon X 2 
representing a resonance arising from interactions between the two defects. In the limit 
X 2 = 0 and (A 2 ,u; 2 ) = (Ai,a;i), these poles appear at 

1 f? 2 A 4 

kp = -2—— *47rA^, (6.16) 

UJi 

In this case, the inverse lifetime has donbled, being the snm of the inverse lifetimes of the 
two oscillators considered separately. It snms for if either of the two oscillators emits a 
photon, the excited state will decay. 

6.2. Band Structure 

In this section we consider a periodic array of defects. With snch a periodic strnctnre, 
we expect the photon spectrnm to be described by a band strnctnre. Rather than compnte 
the band strnctnre by extracting it from Mt, we will be able to compnte it throngh Bloch’s 
theorem. 

We begin with the Hamiltonian describing a three dimensional array of defects in the 
fiber and their interaction with a photon held: 

H = Hfleid + + KS-t + S*t)<Kxi,0,0) + A^4,^(xi,0,0)) . (6.17) 

Here {i,j,k} marks ont the locations of the defects in the array. This array is three 
dimensional. However the problem is still one dimensional as (/>, as before, does not depend 
on the transverse coordinates, y and z. To rednce this Hamiltonian to one dimension, 
we define a collective oscillator representing a snm of all the oscillators in a plane at Xi 
transverse to the direction of the fiber, 

= ( 6 - 18 ) 

^ jk 

where N is the nnmber of oscillators in this plane. has the expected commntator 

IS-,S+] = 1, (6.19) 
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With this commutation relation, the above Hamiltonian can be recast as 


CX) 

= infield + (uJoS+S--^X{S-+S+)(j){ia)+NA^(l)‘^{ia)y (6.20) 


where a is the lattice spacing of the array along the hber. (We do not need to assume the 
lattice spacing in the planes transverse to the hber is also a.) We see both couplings have 
been rescaled, A by y/N and by N. The corresponding Minkowski action is 


S= dtdx 


Stt 




OO 

+ <S(i - li) (-iS."5,S+ - uoS+S" + VN\{S- + S*)4, - NA‘‘4,A 


( 6 . 21 ) 


For a lattice spacing in the transverse planes of 10 ^cm and H ~ 10 ^cm^, we have 
A/wo rsj 10 This is large enough to produce observable results. 

From (|6.21D we can immediately write down the equations of motion governing such 
a periodic array of defects: 


OO r « 

0 = ^ uqNX^ / 


,2 /"“ + (j)n-i{k, (n - l)a) 




2 - F 


1 


- NA'^{4>n{na) + 4>n-i{{n - l)a)) + — [dj:4>n{na) - dj:4>(^ri-i){na)) 
where a is the spacing between the defects. We have written cf) as 

OO 

(l){x) — (l)n{x)9{na — x)9{x — {n — l)a). 


( 6 . 22 ) 


(6.23) 


Each (j)n is free and so satishes (9| — d‘^)(j)n = 0. 

To solve this equation it is necessary to uncouple the various (p^s from one another. 
We do this via Bloch’s theorem. It tells us the solutions to these equations must satisfy 


(j)n{na,t) = e*^“(/)n-i((n - l)a,t). 


This constraint, together with the continuity condition. 


(6.24) 


t') — (pYi — lijXdit')') 


(6.25) 


23 










are enough to specify the band structure. After some algebra, we hnd the consistency of 
(6.20), (6.22), and (6.23), demands 


cos(pa) 


cos(/ca) + 


4nNA^k , 

—5 - TT sin(/ca). 

— Uq 


(6.26) 


Here p is to be understood as the photon momentum and u = \k\ as the photon energy. 



alternatively, underestimated ojqo) in order to r^^Jce the gaps visible. 


A plot illustrating the band structure is given in Figure 1. In this plot are drawn the 
hrst two bands and a portion of the third in the reduced zone scheme. There are three 
notable features to the band structure: 1 ) in the lowest band, the energy, /c, approaches 0 
as /c — 0 ; 2 ) there is a gap between the hrst and the second gap, a result of the resonance 
at /c = Wo; and 3) there is a smaller gap between the second and third bands, more typical 
as it appears as a small perturbation about the edge of the hrst Brillouin zone. 

The existence of a mode at w = /c = 0 results from gauge invariance. Because of the 
simplihcations introduced by wqA^ = A^, the numerator in (6.26) vanishes as /c ^ 0. If the 
term had been omitted from the action, we would have found that long wavelength. 
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low energy photons cannot propagate, as if the photons had an imaginary mass. One, 
however, can ignore the term if, at the same time, one makes the dipole approximation 
and couples the oscillators to the E-held in place of the vector potential, A. This restores 
gauge invariance and so leads to the correct physics. 

The effect of the pole a.t k — uq in (|6.26| ) is to divide what would be the hrst unper¬ 
turbed band k = |/c|, —"Kja < k < ti ja into two bands. The gap between these two bands, 
if small, is given by 


AE = 


47riVA2 


(6.27) 


sin(ct;oti) 

The division into two branches is similar to what one sees in the dispersion relation for 


photons in an covalent crystal (see [0 - we will derive this dispersion relation in the 
next section). However there are two notable differences in Jl^’s treatment from ours. 
There, the defects are coupled to the E-held and not to the vector potential (the dipole 
approximation) and the A^ term is ignored. The second difference is that the pole in is 
shifted away from k = lvq. It seems plausible that this is a result of not treating all the 
oscillators in unison. When we treat an isolated defect in Section 3, we see a similar shift 
in the pole away from uq in T{k) and R{k). 

The last feature of the band structure that bears comment is the gap between the 
second and third bands. This gap, far smaller than the previous one, is equal to 


AE = 




TT^ 


, i2n2 

U/Q Oj 


(6.28) 


We have ignored higher order terms in A^. Similar, but smaller, gaps will appear between 
the higher order bands. Again to lowest order in A^, the gap opening up at coga = rni will 
equal (provided ujq ^ rm ja) 


Ari.E — 


SAtt^A^ 


n 




(6.29) 


Such gaps are what one typically expects to hnd in the case of small perturbations: the 
band structure is only altered around the edges of the Brillouin zones. However there is 
a difference from the standard case: the opening up of the gaps is asymmetric about the 
energy w = rni/a: the two solutions to ( |6.26|) dehning the gap are not w = mi j aE A^E j2 
(symmetric splitting) but w = rni/a and w = rni/a -t- sgn(n^7r^ — ujQa?)AnE. 
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7. Continuum of Defects 


In this section we develop a theory describing a continnnm set of defects living along 
the entire length of the hber. As onr starting point, we take a set of point defects distribnted 
nniformly with linear density p along the hber. The corresponding action is then 


S = j dtdx ^ 

+ - Xk) {-iS-dtS+ - oJoS^S- + \^{S- + S+)(t) - N/\^(j?{xk)) 

k 

(7,1) 

We snppose the rednction to two dimensions, as in (6.17)-(6.20), has already been made. 
To take the continnnm limit of this action, we introdnce the held variables S^{x): 


S^{x) = ~ 

yP k 


(7.2) 


These then obey the expected commntation relations 


[S (x), S~^{x')] = d{x — x'), 


(7,3) 


and the action may be rewritten as 


S 


dtdx 


Stt 




{dx(j))‘^) — iS {x)dtS~^{x) — uoS~^{x)S (x) 


+ ^/I^X{S-{x) + S+{x))(f> - NpA^4>^{x) 


(7.4) 


In doing so, we have smoothed the defect distribntion by nsing the eqnality 5{x — Xk) = 
p, a trick of Dirac’s. As p is typically ~ 10®cm“^, and ~ 10^^, we have A^pA^/wg ~ 10“^. 
This, again, is large enongh to prodnce observable conseqnences. 

From this action we are easily able to derive the exact dispersion relation. As in 
previons sections, we can write down an eqnation of motion involving cf) alone: 


0 = 2 ij 0 qA^N p I duje 


iut . o . 1 


/ >2 , ,2 

(jJ OJn 


+ 2NpA‘^(f){x) + — {df(f) - dl^) . (7.5) 


The mode expansions of cf) take the form 

7 “ 1 / \ 

/ dk _ () 

’—CO 'y'Cj(/c) 7 


(f){x) = 


(7.6) 
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Substituting this into the equation of motion gives us the dispersion relation 


u^{k) = SnNpA^ + + 


2 d>TlN pA'^ojQ 


,2 ■ 


u‘^{k) — Uq 

Like the band structure before, this has two branches. They are described by 


(7.7) 




(i^o + gjrAfpA^ + k^) ± ({k? + SuNpA'^ - + HttN 


(7.8) 


A plot of these two branches is found in the hgure below: 



The modes described by these two branches of the dispersion relation are often referred to 
as polaritons. 

As k ^ oo both branches obtain asymptotic forms, as indicated by the dashed lines. 
The upper branch becomes photonic, u — \k\, as expected. At high energies the photons 
should pass undisturbed through a medium governed by far lower energy scales. The lower 
branch, on the other hand, approaches a constant, u) = loq, in this limit, becoming like an 
optical phonon. As /c —> 0, the upper branch approaches a constant. 


c<;_(_(0) = + d>TiN pA?. 


(7.9) 
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A small gap thus separates these two branches. It is given by /S.E = -^/wq + SnNpA'^ — 
ljq. The lower branch, in contrast, becomes linear as /c —0 and so describes photons 
propagating in a medium of dielectric constant 



(7.10) 


Thus if we know e(0), say experimentally, we can compute the effective coupling strength 


We are also able to compute the frequency dependent dielectric constant that gives 


rise to these dispersion relations. From (|7.8|) and €(u)u‘^ = k'^ we have 



(7.11) 


The dielectric constant, e{u) approaches ±cx) as w —wq. This asymptotic behaviour results 
in the two separate branches of polaritons. As w —> cxo, e{u) approaches 1: as discussed 
previously, at energies far greater than uq the photons cease to see the oscillators. 

8. Impurities in the Presence of a Continuous Medium 

In previous sections, we have treated the defects as if embedded in the vacuum. The 
medium in between the defects has been assumed to have no dielectric properties, i.e. 
its dielectric constant, e{u!), has been set equal to 1. This is not necessarily a bad ap¬ 
proximation. The material out of which a hber is constructed typically has a dielectric 
constant with little frequency dependence, i.e. e{u) — eg, for a wide range of frequencies, 
hujQ < lOeV. So by a judicious insertion of cq in the equations from previous sections, we 
can take into account the medium of propagation. However if we wish we can do better. As 
we have developed here a theory of a continuous medium, it is possible to solve the problem 
of a defect embedded in a medium with a frequency dependent dielectric constant. We 
will consider the defect to arise from atomic polarization. Alternatively we could consider 
the defect to be ionic in nature and suppose the natural frequency of the medium to be 
far below the energy scale where the ionic defect will also experience atomic polarization. 
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The action describing snch a defect is a sum of medium’s action and the defect’s 
action: 


S — 5'nied + Sd 


1 


— ({dt4)Y - - iS {x)dtS+{x) - LJmS^{x)S (x) 


= J dtdx 

+ \/NmPm>^m{S~{x) + {x))(j) - 


( 8 . 1 ) 


- j dt [iS-dtS+ + ujdS+S- - Xd{S- + S+)(j) + A2(/,2] . 

We have labeled the defect and medium couplings with d and m subscripts respectively. 
The only difference in solving this theory from previous ones is in the nature of the mode 
expansions. Here the mode expansions are done in terms of polaritons: 


j=+,-' 


f*00 


1 


/ dk- ^ _ 

-oo \/wj{k) 




( 8 . 2 ) 


The sum reflects that the dispersion relation has two branches. So the scattering matrices 
for the polaritons off the defect are similar to those of (|3.12|) with u!{k) = \k\ replaced by 
u!{k) = u±{k) (where u!± are given from ( [T^ ) by replacing ojq with Np with NmPm, 
and A with A^,). Thus we have 


T±{k) 

R±{k) 


ik{i^±{k) - ul) 


ik{oj‘i{k) - uf) - AnAlulik) ’ 
_ 47rA^u;|(fc) _ 

ik{ujl{k) - uj) - 47rA2w^(/c)' 


(8.3) 


We now need to index T and R to indicate which branch the polariton is on. 

Similarly we can imagine an array of defects imbedded in the continuous medium. 
The band structure is then derived from ( |6.26|) by replacing the free dispersion relation 
with ij±: 

47rA^.AT,)2 (k'] 

(8.4) 


47rAdA^Ci;i(/c) 

cos(p±a) = cos(/ca) + ^ ^ 2 n \ - ^sin{ka), 


k{u^Uk)-^^d) 

where Nd is the number of defects in a section transverse to the hber. Here p± is the 
crystal momentum for polaritons on the upper/lower branch of the dispersion relation. 
For a given p±, a set of /c’s, {/cj}, satisfy this equation. The energies corresponding to this 
p± are then {u!±{ki)}. A plot of the band structure is given in Figure 3 below. 
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Figure 3: A plot of the band strijcture for a periodic of defects imbedded in a 
material governed by the dispersion relation in (7.7). For the defects we have again 



taken = 2(^47rA^dC<;^A^a^ =^. For the parafheters of the medium, we have 

taken Uma = 6 and 327rAmPmA^c<;^a^ = 100. xB^gaps that are too small to be 
seen have been explicitly marked out. 




4 


The gap between uj+{k) and U-{k) forces the bands corresponding to these two 
branches to be separated by this same gap. The bands of U-{k) are found below Um 
while those of uj+{k) are bounded below by The lower set of bands 

away from LVm are similar to those found in Section 6. For small k and w, the spectrum 
is linear; there is a large gap between the hrst and second bands; and there are smaller 
gaps between the higher order bands. The similarity results from uj-{k) being photonic for 
small k. Close to ujm the bands of uj-{k) accumulate. They are of an inhnite number, each 
separated by vanishingly small energies. The existence of such bands owes itself to the 
asymptotic form of u-{k): as /c —> cx), U- approaches a constant. The upper set of bands 
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are photonic away from Um (provided Um > ^d) as rc+(/c) is photonic for large energies. 
Only near Um do the bands of c<;+(/c) deviate from linearity, a resnlt of Ci;_|_(0) 7 ^ 0. 

The gap between the hrst and second band of U-{k) is described as in (6.27): 



where e{ud) is given by (7.11). The gaps opening np at p = rnr for both sets of bands are 
given by 

AnE± = \'y±{n)\dkU±{n7r/a)-, 


u+{k) 

dkUJ+{k) 


_ 8iVdA^w|(n7r/a) 
n{ul{mT/a)-uj)' 
k+0{A‘^); u)_{k) = u)m + 0{A‘^)-, 


1+0(A^); dku-{k) 


SttA^A^ Pm^mk 

(/c2-w^)2 


( 8 . 6 ) 


We have again only kept terms to lowest order in O(A^) and O(A^). These gaps, as 
with their connterparts in Section 6, are not symmetric abont u±{n7T). Rather the gap is 
dehned by w = c<;±(n7r/a) and u = u)±(nTi/a) + sgn(a;±(n7r/a)^ — uj'^)AnE±. 
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